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ABSTRACT 

The complete set of Higgs-boson two-loop contributions to electric dipole moments of the 
electron and neutron is calculated in the minimal supersymmetric standard model. The 
electric dipole moments are induced by CP- violating trilinear couplings of the 'CP-odd' 
and charged Higgs bosons to the scalar top and bottom quarks. Numerical estimates of 
the individual two-loop contributions to electric dipole moments are given. 
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Supersymmetric (SUSY) theories, including the minimal supersymmetric standard 
model (MSSM), face the difficulty of explaining naturally the apparent absence of electric 
dipole moments (EDMs) of the neutron and electron ||T|, ^, |^. Several suggestions 
have been made to suppress the SUSY contributions to electron and neutron EDMs, at a 
level just below their present experimental 2cr upper bounds: \df.\ < 0.5 x lO^^^e cm and 
\dn\ < 1-12 X 10~^^e cm [||. Apart from the obvious choice of suppressing the new CP- 
violating phases of the theory to the 10"'^ level [|l|, 0, a more phenomenologically appealing 
possibility is to make the ffist two generations of scalar fermions as heavy as few TeV, 
but keep the soft-breaking mass parameters of the third generation relatively small, e.g. 
0.5-0.7 TeV 0. An interesting alternative is to arrange for partial cancellations among the 
different EDM contributions 0], within the framework of superstring-derived models [0]. 

Here, we shall focus our interest on studying additional two-loop contributions to 
electron and neutron EDMs in a SUSY scenario, in which the ffist two generations are 
rather heavy, e.g. of order few TeV, whereas the third generation is relatively light below 
the TeV scale [Q. In such a theoretical framework, the leading effects on the EDMs arise 
from CP-violating trilinear interactions related to the scalar top and bottom quarks through 
the three-gluon operator and through the coupling of the 'CP-odd' Higgs boson, a, to 
the gauge bosons , as shown in Fig. |^. Despite their similarity to the graphs due to Barr 
and Zee 0, our reasoning of considering the two-loop contributions of Fig. ^ is completely 
different from [^J, as the third-generation scalar quarks can have a significant impact by 
themselves on the electron and neutron EDMs, independently of the chirality-unsuppressed 
one-loop contributions. In our analysis, we shall assume rather heavy gluino masses m^, 
e.g. rrig > 0.5 TeV, such that their effect through the three-gluon operator, which scales as 
1/m? 0, will be much smaller than the present experimental upper bound [Q. 

In this paper, we shall complement the analysis of Ref . ^ , and consider the complete 
set of two-loop EDM graphs, shown in Fig. |l|, including those due to CP-violating aA^Zx 
and if"'"y4^W3(7|^ couplings. At this point, it is worth stressing that the EDM constraints 
we shall study here will have important consequences on Higgs-sector CP violation within 
the MSSM found recently and on related phenomena in 5-meson decays, dark-matter 
searches and collider experiments [0 . 



Our starting point is the scalar top and bottom mass matrices, which may conve- 



* Analogous EDM contributions induced by charged-Higgs-boson two- loop graphs were studied in [|IC 
within the context of three-Higgs-doublet models. 
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niently be expressed, in the weak basis (g^, qn), as follows: 



Ml 



+ ml + cos 2/?M| (Tf - Q, sin^ 9^) m,(A* - fxRg) 

mg{Aq - ij,*Rg) + ml + cos 2/5M| Qq sin^ 

(1) 

with q = t,b, Qt = 2/3, Q,, = -1/3, = -T^" = 1/2, Rb = tan(3 = V2/V1, Rt = cot (3, 
and sin6'^ = (1 — M^/M|)^/^. Moreover, Mq and are soft-SUSY-breaking masses 
for the left-handed and right-handed scalar top and bottom quarks. The matrix Ad^ may 
be diagonalized through a unitary rotation, which relates the weak {q^, qji) to the mass 
eigenstates (gi, ^2): 



Ql \ _ ( ^ W cos6'g sinOq \ I qi 
Qr I \ e'^" ) \ -sin^g cos^^ I \ 



(2) 



where 5q = arg(Ag — Rqjj*) and 6q are mixing angles determined by 

cos On 



mq\Aq- RqH* 



ml\Aq-Rq^i*\^ + [{Ml)LL-Mlf 



■ n \{MI)ll-MI\ 

sm.6q = ^ = . (3) 

"\Aq-Rq^^*\^ + [{MDlL-MI]^ 



m„ 



In Eq. @, the quantity {M.'^^ll is the {ll}-matrix element of the scalar quark mass matrix 
M.^. In addition, the mass eigenvalues of M.^ are given by 



Kim) = ^ {m2 + M^' + 2mJ + T,^ cos 2/3M| 



+ 



- M2 + cos2/?M|(T| - 2Qq sin^ e^)f + 4ml\A*q - fiRq\^ |. (4) 



As usual, we consider the convention in which arg/x and argA^ are the only physical 
SUSY CP-violating angles in the MSSM. Then, CP violation originates from the interaction 
Lagrangian 

Ccp = J2 {q*ReT'''^*'^^q,) + {h^ i* tlmV'''''*'^ b, + H.c.) . (5) 

q=i,b *.i = 1.2 «J = 1,2 

The real and imaginary parts of the couplings aq*q and H~^t*bj, which are denoted by r°-^i^J 
and r^^^i^^ , respectively, are given by 

' ' -cot2^„ -1 ' ^ ^ 
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/ cosgt c _|_ cosgj, c co&Qt C _|_ sing;, c 

jj^-p-ff+t*bj _ J__ I cosOt^^ ~^ cos 6»t sinOt '^^ "T cose*'^' 

I sin 9t c cos Of, p sin Qt C sin Q\y /• 

\ cos6»b sine* ^* sinSft sin 64 ^' 



^ -\/2mbm4 sin(5{, — 5t) / sin6'fsin6'fc — sin cos i 
sin/?cos/3f2 I - cos 6't sin 6^6 cos6'tcos^b 



where v = \jv\\ v\ = 2gyj/Mw, and and are the CP-violating quantities 

sin26'gmglm(/ie*''9) 
sm p cos p f 

We shall first calculate the couplings aA^j_{k)A\{q), agfj_{k)gx{q), aA^{k)Zx{q) and 
H~' Afj^{k)W^ {q) , which are induced by t- and 6-mediated one-loop graphs shown in Fig. 
|I|. We adopt the Feynman-'t Hooft gauge, and neglect Feynman diagrams that lead to 
suppressed EDM contributions proportional to m'j. With the convention that the momenta 
q and k fiow into the vertex, the analytic result of the one-loop couplings is found to have 
the gauge-invariant form 

tTl^{k,q) = [{q-k)g^x - q^^h] , (9) 

where the superscript V denotes the gauge boson propagating in the loop, and 

X / dx — ^ TT-r^ TT—, r , (11) 

io Mix + Ml{l - x) ~ q^x{l - x) ' ^ ' 

aW(,2^ _ Nc9w(^ ^ frH+ifbA* r.tb x{l-x)[QtX + Qbjl-x)] 

^ " Sv^TT^ ,,jit,2 ^ ^ ' Jo Mfx + MfXl -x)- q'x{l - x) '^''^ 

The corresponding one-loop form factor A^ may be obtained by replacing the colour factor 
AT^ = 3 by 1/2, and by 1 in Eq. (|10]). In Eqs. ([ll]) and (|T|), {q = t,b) and K^^ 
are (2 x 2)-non-unitary matrices describing the mixing in the Zq*qj and W~^t*bj sectors, 
respectively: 



^ ( 2T,^cos2e,-2Q,sin2^^ T,«sin2e, . 

' Tfsin2^g 2T,« sin^ - 2Qg sin^ ^ ^ 



j^tb _ I cos 9t cos 9b cosOtsin 9b , ^^^^ 
sin 9t cos 9b sin 9t sin 9b 



4 




Figure 1: Higgs-boson two- loop contributions to EDM and CEDM of a fermion in the 
MSSM in the Feynman-'t Hooft gauge (mirror graphs are not shown); /' represents the 
conjugate fermion of / under T/. 

At this stage, we should remark that there are also contributions originating from chargino 
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loops. However, these contributions are proportional to arg/i, on which strict constraints 
exist from one-loop graphs that contribute to the electron EDM, and may therefore be 
neglected. Finally, we should notice that the one-loop AW ^-coupling receives its gauge- 
invariant form of Eq. @ in the Feynman-'t Hooft gauge, after including the H~W~^ and 
H^G^ wave functions and the respective tadpole contribution related to the H^G^ tran- 
sition, as shown in Figs. |l|(e)-(g). 

It is now straightforward to compute the individual contributions to the EDM of a 
fermion that come from quantum corrections involving 7, Z and bosons in the loop. 
These individual EDM contributions shown in Fig. |l] may conveniently be cast into the 
form: 



dj 



w 



N^a tan/3 m/ 



q=t,b 



M? 

gi 
M2 

a 

tan/5 rrif 



64 cos^ 6,,, tt'^ 



V V - BeT"^'^' K'^ O g(— ^ 

l=t,b ij=l,2 ^ a a 

tan/5 rrif 



Ml 
2 Ml 



1 



128V27r3 



QtG 



E 



'■J 



Ml, 



M? 



Ml. 



Ml 
Ml. 



mIl 

Mh 



M2 



M, 



Ml J 



Ml 

Ml. 



M? 
Ml. 



(15) 



(16) 



(17) 



where a = e^/(47r) and = g'^/{ATT) are the electromagnetic and weak fine structure 
constants, respectively, and F{z) and G{a, b, c) are two-loop functions given by 



Fiz) 



G{a, b, c) 



+ 



dx 



x(l 



In 

axil 



xil — x) 



z 

x) In a 



a — l)[ax(l — x) —bx — c(l — x)] 
x(l — x)\px + c(l — x)] 



[ax{l — x) — bx — c(l — x)] [x{l — x) ~bx — c(l — x) 



In 



/bx + c(l — x) 
\ xil — x) 



(19) 



Note that 2G(0,a,a) = -F(a). Among the EDM terms given by Eqs. ([T|)-(|T3), the 
fermion EDM induced by the photon-exchange graphs {df/ey represents the dominant 
contribution For completeness, we give the two-loop contribution to the CEDM of a 
coloured fermion 



(ij \ s as tan /? mj 



647r3 



M2 

a 



q=t,b 



M? 

qi 
M2 J 

a ' 



(20) 
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We can now estimate the neutron EDM dn induced by du and dd in the valence quark 
model, by including QCD renormalization effects H, i.e. 



dn 

e 



dd 

e / A 



d^ 



(21) 



The u- and d-quark masses occurring in Eq. (|T^) are running masses evaluated at the low- 



energy hadronic scale A. In Eq. (21), we have assumed that the renormalization-group 
running factors of the strong coupling constant from mf, to A is almost of order 1. To 
be specific, we consider the values: m„(A) = 7 MeV, md{A) = 10 MeV, as{Mz) = 0.12, 
and 5's(A)/(47r) = I/Vg 0. Likewise, the light-quark CEDMs d^ and d^ lead to a neutron 
EDM 

gs{Mz) 



d^ 



9s{A) 



9 



+ 



2 /d' 



9 



c 

"u 

e 



(22) 



where in turn the strong coupling constant Qg and the u- and ci-quark masses in d^ and d^ 
are calculated at the scale A. 

In the following, we shall give a more quantitative discussion of the individual two- 
loop contributions to electron and neutron EDMs for Ma = 150 and 300 GeV. At the tree 
level, the charged-Higgs-boson mass Mh+ is related to the would-be CP-odd mass Ma by 



M'a + Ml 



(23) 



Even though this very last relation receives appreciable radiative corrections in the MSSM 
with Higgs-sector CP violation |T^, we shall still make use of Eq. (PBj), as required by a 
consistent expansion in perturbation theory. As we have explicitly demonstrated in [^], the 
EDMs crucially depend on /i and tan [3 through ^t,b in Eq- (H) and through the couplings of 
a and to electron and d quark. For the purpose of illustration, we therefore plot in Figs. 
and the numerical predictions for electron and neutron EDMs as functions of tan (3 and 
/i, respectively. Specifically, we consider the following values for the SUSY parameters: 



Fig. 2 : Mo = Mq = M^ = Md = 0.6 TeV, A = \At\ 



I A 



A = ^ = I TeV, argA = 90° 
Fig. |: Mo = Mq = M^ = Md = 0.6 TeV, A = |A| = |A| 
tan/? = 20, arg^ = 90° 



1 TeV, 



1 TeV, 



(24) 



In addition, we assume that the yU-parameter is real. In agreement with [Q, we find that 
the dominant EDM effects originate from the aA^AA-coupling in (ig, {de)"' ■, and from the 
ag^g\-co\x\Am.g in dn, (c^n)*"- In Figs. ||(a) and|^(a), we see that the charged-Higgs-boson 
two-loop contribution to EDM, ((ig)^, is smaller than {d^.)'^ by a factor 8. On the other 
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hand, Figs. 0(b) and |^(b) show that the corresponding EDM contributions, (c?„)^, (dn)^, 
and (dn)^ , are all of comparable size. They are roughly one order of magnitude smaller 
than (dn)'", and of opposite sign. 

In conclusion, we have shown that EDM constraints on the CP-violating parameters 
related to the sectors of scalar top and bottom quarks can be significant for tan/? > 10, 
and fi, At^h - 0.5 TeV. In particular, we have studied additional two-loop contributions 
mediated by W- and Z-boson interactions, which are found to be sub-dominant but non- 
negligible, and are therefore expected to play an important role in future phenomenological 
analyses of Higgs-sector CP violation in the MSSM. 



Note added 

While revising the paper, I became aware of Ref. which addresses the same topic. 



After the final revisions, the results obtained by the two groups agree both analytically 
and numerically. The author wishes to thank Darwin Chang and Wai-Yee Keung for 
discussions. 
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Figure 2: Numerical estimates of the individual two-loop EDM contributions as a function 
of tan/?: (a) (de)'^ (solid line), (de)^ (dashed line); (b) —{dn)'^ (solid line), (dn)'^ (dashed 
line), (dn)^ (dotted line), ((i„)^ (dash-dotted hne). Lines of the same type from the upper 
to the lower one correspond to = 150 and 300 GeV, respectively. 
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Figure 3: Numerical estimates of the individual two-loop EDM contributions as a function 
of /x: (a) {de)'^ (solid line), {de)^ (dashed line); (b) —{dn)^ (solid line), {dn)'^ (dashed 
line), {dnf^ (dotted line), ((i„)^ (dash-dotted hne). Lines of the same type from the upper 
to the lower one correspond to = 150 and 300 GeV, respectively. 
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